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Abstract 

Rapidly rotating Myers-Perry black holes in d > 6 dimensions were conjectured to be unsta- 
ble by Emparan and Myers. In a previous publication, we found numerically the onset of the ax- 
isymmetric ultraspinning instability in the singly-spinning Myers-Perry black hole in d = 7, 8, 9. 
This threshold signals also a bifurcation to new branches of axisymmetric solutions with pinched 
horizons that arc conjectured to connect to the black ring, black Saturn and other families in the 
phase diagram of stationary solutions. We firmly establish that this instability is also present in 
d = 6 and in d = 10, 11. The boundary conditions of the perturbations are discussed in detail 
for the first time and we prove that they preserve the angular velocity and temperature of the 
original Myers-Perry black hole. This property is fundamental to establish a thermodynamic 
necessary condition for the existence of this instability in general rotating backgrounds. We 
also prove a previous claim that the ultraspinning modes cannot be pure gauge modes. Finally 
we find new ultraspinning Gregory-Lafiammc instabilities of rotating black strings and branes 
that appear exactly at the critical rotation predicted by the aforementioned thermodynamic 
criterium. The latter is a refinement of the Gubser-Mitra conjecture. 
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1 Introduction 

The turn of the century witnessed a renewed interest in Einstein's gravity. The research focus 
extends nowadays to numerical simulations of the full time evolution (inspiral, merger and ring- 
down phases) of black hole binary systems [1], numerical simulations of high energy collisions 
of black holes [2], and the study of higher-dimensional black holes [3], mainly in vacuum and 
in asymptotically anti-de Sitter (AdS) spacetimes. The former two endeavours are of the utmost 
interest for gravitational wave experiments. On the other hand, the latter programme was triggered 
by the quest for a microscopic description of black holes in a theory of quantum gravity (string 
theory), by the emergence of TeV-scale gravity scenarios relevant for the LHC, and by the realization 
that black holes describe thermal phases of gauge theories in holographic gauge/gravity dualities. 
The properties of higher-dimensional black holes, namely their stability, will be the main focus of 
the present work, and we review its motivations in more detail next. 

Black hole thermodynamics and the thermal spectrum of Hawking radiation strongly suggests 
the existence of a statistical mechanical description of black holes in terms of some underlying 
microscopic degrees of freedom. This microscopic description necessarily requires quantum gravity. 
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One of the most compelling candidates for such a theory is string theory, which inevitably requires 
gravity in extra dimensions. Furthermore, the observation that the black hole entropy is propor- 
tional to the horizon area strongly suggests that the quantum degrees of freedom of the black hole 
are distributed over a surface [3j . This led to the formulation of the holographic principle according 
to which quantum gravity in a given volume should have a description in terms of a quantum field 
theory on its boundary. The last decade witnessed concrete realizations of this holographic idea, 
namely the AdS/CFT duality and more generically the gauge/gravity correspondence The ap- 
plication of these ideas to certain classes of higher-dimensional black holes has already provided a 
statistical description of some of their thermodynamic properties. Successful examples include: i) a 
statistical counting of the Bekenstein-Hawking entropy [6] , ii) a microscopic description of Hawking 
emission in near-extremal black holes [7] , iii) a microscopic description of superradiant emission [8] , 
iv) the map between the Hawking-Page phase transition in asymptotically AdS black holes and the 
confinement/deconfinement phase transition in gauge theories that share some common properties 
with QCD [9], and v) the identification of the quasinormal modes of oscillation of a black hole with 
the thermalization frequencies of the perturbed holographic field theory jlOj , to mention but a few. 
Following these early successes, the gauge/gravity correspondence has matured and is nowadays a 
field of research on its own. Mainly due to the fact that this is a weak/strong coupling duality, it 
is recognized as a powerful technical tool to understand not only gravity but also strongly coupled 
gauge theories. In recent years this correspondence has been extended to other systems of interest, 
namely, the fluid/gravity correspondence I12j and more recently the condensed matter /gravity 
correspondence jl3| . These developments indicate that gravitational physics is becoming a valuable 
technical tool to understand several other branches of physics. Given the holographic nature of 
these correspondences, and since we are often interested in gauge theories in four dimensions, it is 
therefore important to develop our understanding of gravity in five and higher dimensions. 

The LHC will soon be operating at the TeV energy frontier. Although it was designed to find 
the Higgs boson and study particle physics beyond the Standard Model, there is the possibility 
that the LHC discovers extra dimensions. Indeed, in recent years various TeV-scale gravity or 
braneworld scenarios have been proposed according to which extra dimensions might be detected 
at the energy scales probed by the LHC (see |14j for a review). The motivation for these proposals 
is the solution of the hierarchy problem, i.e., the huge difference between the Planck scale and 
the electroweak scale. By proposing the existence of sufficiently "large" extra dimensions, the 
fundamental (d-dimensional) Planck mass can be of the order of the electroweak energy (~ 1 TeV) . 
The most fascinating outcome of these scenarios is the possibility of producing microscopic black 
holes at the LHC. This could allow for the experimental determination of the fundamental scale of 
gravity, the number of extra dimensions, and the decay product of higher-dimensional black holes, 
including the first evidence of Hawking evaporation [T3]. This fascinating possiblity provides an 
extra motivation to undergo the study of higher-dimensional black holes. 

Stimulated by the advances in higher-dimensional gravitational physics, there is an intensive 
ongoing program to extend analytical tools and numerical relativity to higher dimensions. The 
Newman- Penrose and Geroch-Held- Penrose formalisms, and the resulting Petrov classification, 
which are important to classify solutions and study their stability, have been extended to d > 4 [T5] . 
The separability of some class of perturbation equations and the existence of completely integrable 
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geodesic equations in some higher dimensional black holes is possible due to the fact that there 
are hidden and explicit symmetries associated with the principal conformal Killing- Yano tensor 
|16j . The study of higher-dimensional gravity is also motivating the development of completely 
new analytical tools, like the blackfold approach [n\ I18j. to construct perturbative solutions based 
on the fact that for higher-dimensional black holes the horizons can have widely separated length 
scales. Simultaneously, the first non-linear numerical codes in higher dimensions are starting to 
be developed mainly to study the time evolution of instabilities in black holes |19[ \2T\ [20] . This is 
expected to be an area of active research in the near future. In addition, the challenges of higher- 
dimensional black hole physics have motivated the introduction of new numerical tools, like the 
spectral methods [22l[23l[2i]. 

Last but not least, another motivation to study higher-dimensional black holes is the under- 
standing of these objects per se. Indeed, to uncover the mathematical structure of Einstein's gravity 
and its solutions one should treat the number of spacetime dimensions d as a free parameter in the 
theory. We should be able to distinguish between the universal properties of the theory and the 
dimension dependent ones. 

It is known since the 70's that 4(i black holes are the simplest gravitational objects in Nature. 
Hawking's black hole topology theorem states that a 4d black hole must have an event horizon with 
spherical topology [25]. Together with Hawking's rigidity theorem |25| this led to the proof of the 
4d uniqueness theorems ^Ei- As a result, a 4(i vacuum black hole is fully specified by its conserved 
charges: mass M and angular momentum J. Therefore, there is a unique black hole in 4d Einstein 
gravity — the Kerr solution — which has an upper bound on its angular momentum, J < GM"^ 
{G being Newton's constant). However, are these universal properties of gravity? Or are some of 
them dimension dependent? 

The last decade has provided a cascade of fascinating answers to these questions (see [3] for 
a review). The higher-dimensional arena is populated by new solutions besides the Myers-Perry 
(MP) black hole with S'^^'^ topology [27] — the (non-trivial) counterpart of the Ad Kerr solution. 
The simplest black object in d dimensions is a black string with horizon topology S"^^^ x R or a 
black brane with S''^~^~" x M", constructed by adding trivial flat direction(s) to the metric of the 
Schwarzchild geometry. The most surprising analytical solution is however the bd black ring found 
by Emparan and Reall [28]. This solution is asymptotically flat, has horizon topology S'^ x 
and in some regions of the parameter space it can carry the same conserved charges as the singly- 
spinning Myers-Perry black hole. Therefore, this solution explicitly demonstrates that the topology 
and uniqueness theorems do not generalize straightforwardly to higher dimensions when rotation 
is considered!^ Indeed, the generalisation of the black hole topology theorem imposes much weaker 
restrictions on the horizon topology [31]. In addition, the 5d uniqueness theorems, restricted to 
geometries with two U{1) Killing isometrics, explicitly show that in order to uniquely specify a 
black hole solution, one has to fix the conserved charges and other parameters not related to them 
[32j (see also [33] for some earlier work). Using the complete integrability of the five-dimensional 
Einstein's vacuum equations restricted to solutions with M x f/(l)^ Killing isometrics, it has been 
also possible to construct explicitly the black rings with rotation only along the S'^ (which possess 

^The situation is considerably different for static solutions. Indeed, it is proven that the Schwarzschild-Tangherlini 
black hole is the unique static solution [2S], and moreover this black hole is stable at the linear mode level |30j . 
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conical singularities) (Mjj^tlie regular doubly-spinning black ring ^36j, and regular asymptotically 
flat multi-black hole objects like black Saturns [37j, di-rings [38] and bi-cycling rings [39]. These 
solutions provide further examples of the richness of higher d gravity. In d > 5, the integrability 
methods are not available, but the perturbative blackfold approach of [T7j has provided evidence 
for the existence of black objects with other (non-spherical) topologies [l0l[18]. All these solutions 
have multiple U{1) rotational symmetries and need to be organized in a classification scheme. To 
provide such a classification we must address the question of to what extent the rigidity theorem 
extends to higher dimensions. In higher dimensions, stationarity only implies the existence of one 
rigid rotational symmetry |41j. However all higher-dimensional solutions known exactly have more 
than one such symmetry. A fundamental question is whether there are stationary black holes with 
less symmetry than those solutions. This possibility was first raised in |42j . where the existence 
of stationary black holes with a single U{1) was conjectured. We shall return to this fundamental 
issue below. 

Having these explicit solutions, the next natural question is whether they are dynamically stable 
against linear perturbations. Start with the Kerr solution in vacuum. Smarr observed that, as the 
rotation increases, the horizon Gaussian curvature at the poles starts positive, becomes zero at a 
critical rotation, and then goes negative before the Kerr bound J = GM'^ is reached [H]. He further 
noticed that a similar behaviour occurs in rotating fiuid droplets held by surface tension, where 
this behaviour signals an instability. However, soon after Smarr made this remark, Teukolsky 
used the Newman-Penrose formalism to find the master equation that governs the gravitational 
perturbations of the Kerr black hole, which was used by Whiting to show the mode-by-mode 
stability of the Kerr black hole [45j. 

What about MP black holes? For d > 6, the rotation of a singly-spinning MP can grow 
unbounded, as already noticed in [27]. This led Emparan and Myers to propose that, for sufficiently 
high rotation, these black holes should become unstable against what was called the ultraspinning 
limit [46]. The natural way to check the presence of this instability would be to apply a linear 
perturbation analysis using the Newman-Penrose formalism, this time for d > 6. Unfortunately, 
although this formalism and the associated Geroch-Held-Penrose formalism have been extended 
to higher dimensions, it is extremely difficult if not impossible to manipulate the equations to get 
decoupled master equations for the perturbations [l5]. So, in practice, we cannot use it yet to 
treat the perturbation problem analytically]^ Given these limitations, Emparan- Myers provided 
solid heuristic arguments, that we highlight next, to conjecture the existence of an instability for 
a sufficiently large rotation [46j. Take a MP black hole rotating along a single plane for simplicity. 
In the ultraspinning regime, a ^ r+ (where a = ^^jj, and r_|_ is the horizon radius), the black 
hole horizon fiattens out along the plane of rotation and its shape can be approximated by that 
of black disk of radius a and thickness r+ |46| . In fact, one can take a precise limit in which the 
MP black hole metric near the rotation axis reduces to the metric of a static black membrane with 
horizon topology S'^~^ x M?, and where the worldvolume directions of the membrane correspond to 

^Ref. |35j was able to construct the same solution in "ring-like" coordinates without using integrability methods. 

''Some subsectors of the perturbations of some classes of MP black holes have been decoupled [491 150| . but none 
of them shows signs of any instability and indeed they do not contain the kind of perturbations relevant for the 
ultraspinning instability. 
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the directions along the original rotation plane |46j . This is a far-reaching observation since black 
membranes are afflicted by the Gregory-Laflamme (GL) instability [17j. This led the authors of |l6] 
to conjecture that ultraspinning MP black holes should be unstable under a Gregory-Laflamme- 
type of instability (see subsection 2.1 for more details). According to the arguments of [l6], this 
instability should become active when a ~ r-(_, and it provides an effective dynamical bound on the 
angular momentum. As often occurs in physical systems, and in particular in the context of black 
branes [48j, the stationary threshold of the ultraspinning instability is expected to signal also a 
bifurcation point in a phase diagram of solutions into a new branch of pinched black objects. New 
thresholds are expected at higher rotations, from which additional branches of pinched solutions 
should bifurcate. These conjectured stationary solutions preserve the original MP symmetries, 
but have spherical horizons distorted by ripples along the polar direction, and are conjectured to 
connect to the black ring and black Saturn families in a phase diagram of stationary solutions [40]. 

The main purpose of the present study is to firmly establish the properties of the threshold of this 
ultraspinning instability. A practical use of the Newman- Penrose formalism in higher dimensions 
is not yet available so we will resort to a numerical approach. In a previous publication, we have 
already addressed this problem in d = 7,8,9 [23j. The present manuscript, however, extends the 
former study in several directions. We present for the first time the properties of the instability 
in the six-dimensional case. This is an important case since d = 6 is the lowest dimension where 
the instability is present, and it was not addressed in |22i due to numerical difficulties. Moreover, 
we discuss the instability in d = 10 and d = 11 dimensions that are relevant for string theory. We 
also address the d = 5 case, in which the black hole does not have an ultraspinning regime (in 
the sense that defined in subsection 2.2). Consequently, no instability that preserves the spatial 
isometrics of the background can be present, as we confirm numerically. At a more technical level, 
we will provide the first detailed discussion of the boundary conditions of the problem. While doing 
so we will prove that the ultraspinning modes preserve the temperature and angular velocity of 
the background, an aspect that was not discussed in [23]. This is a keypoint to formally establish 
the thermodynamic criterium (which is a necessary condition, not a sufficient one) for the critical 
rotation above which black holes can be afflicted by the ultraspinning instability. It can be seen 
as a refinement of the Gubser-Mitra conjecture |51l I52j . This thermodynamic criterium was first 
proposed in (23], and further developed in [21]. We will also prove the claim of [23] that the 
ultraspinning threshold modes that we find cannot be pure gauge modes. 

The ultraspinning instability is not unique to singly-spinning MP black holes. Indeed, [23] found 
that it is also present in the case of MP black holes with equal angular momenta along the [{d—1) /2j 
rotation planes that are allowed in d dimensions (here [J stands for the smallest integer part). This 
occurs in spite of the fact that for these black holes, contrary to the singly-spinning case, the 
angular momenta have an upper bound. The keypoint is that the aforementioned thermodynamic 
criterium (subsection 2.2) still allows for an ultraspinning regime. In the equal spins case and 
in odd spacetime dimensions, the geometry is codimension-1, i.e., it only depends on the radial 
coordinate. This simplifies considerably the analysis of the perturbed Einstein equations since they 
reduce to a coupled system of ODEs. For this reason, the detailed study of the time dependence 
of the instability is much simpler, and |[24j found, as expected, that the modes responsible for the 
ultraspinning instability grow exponentially with time. This fact provides strong confidence to 
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expect that a similar situation occurs in the singly-spinning case. 

The analysis of the ultraspinning instability in the equal spins case [24] gave the first evidence 
for the existence of a new family of topologically spherical black hole solutions with only a single 
rotational symmetry, i.e., solutions that saturate Hawking's rigidity theorem in higher dimensions 
|41j . There is a previous example of a black object with a single rotational symmetry but different 
horizon topology [18]: the "helical" black rings that can be perturbatively constructed using the 
blackfold approach |17| . These developments confirm the conjecture put forward in [32] and can 
fairly be seen as one of the most important legacies of the ultrapinning instability studies and of 
the blackfold approach. More recently [13], near horizon geometries of extremal black holes have 
been found with a single rotational symmetry. 

The plan of the paper is the following. In section [2| we will review the ultraspinning instability 
conjecture as formulated in [36], together with the refined thermodynamic criterium of |23t I24j . and 
we will formulate the problem of finding the modes responsible for the ultraspinning instability as 
an eigenvalue problem. In section [3j we will discuss the boundary conditions appropriate to the 
problem at hand and we will prove that the threshold modes searched numerically cannot be pure 
gauge. In section |4j we will present our results and finally we will close with a Discussion. An 
Appendix contains the technical details of the horizon embeddings. 



2 The ultraspinning instability 



In subsection|2.1[ we will review the work of [46J where the conjecture of the ultraspinning instability 



of singly-spinning MP black holes was first formulated. In subsection 2.2, we will review the 
thermodynamic arguments that led the authors of [23\ [23] to extend the ultraspinning conjecture 
to other classes of black holes which need not admit an arbitrarily large angular momentum (per 
unit mass). Finally, in subsection 2.3, we will formulate our perturbation problem as an eigenvalue 
problem. 



2.1 Motivations 

The metric of the d-dimensional asymptotically fiat MP black hole spinning in a single plane is 
given by [27] 

ds' = - {dt + a sin^ 9dc^f + [(r^ + a')d^-adt\' + dr' + S(r, 6) dO' 

h[r,U) h[r,U) '- -' A[r) 

+ cos^ 9 , 

(2.1) 

where dQ,'^^_^^ is the line element of a unit {d — 4)— sphere and 

j:{r,e) = r^ + a^cos'^9, A(r) = + - ^ . (2.2) 

This solution of the Einstein vacuum equations is characterised by two parameters, namely the 
mass-radius and the rotation parameter a, 

rf„3 _ WttGM _d-2 J 

"■"^ ~ {d- 2)Ad-2 ' ""~M' ^ ^ ^ 
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where Ad-2 = 2'ir^^~^'^/'^ /T[{d — l)/2] is the volume of a unit-radius {d — 2)-sphere, G denotes 
Newton's constant and M and J are the ADM mass and angular momentum, respectively. The 
(outer) event horizon lies at the largest real root r = of A(r) = 0, that is, 

m 

d-b 



rl + a'---^ = Q. (2.4) 



For d = 4 a regular horizon exists for all values of a up to the Kerr bound, a = rml2 (a = GM), 
which corresponds to an extremal (i.e., zero temperature) black hole with finite size horizon area. In 
d = 5 the situation is somewhat similar: a regular solution exists for all values of a strictly smaller 
than r^j^ The solution with a = corresponds to a naked singularity. On the other hand, for 
d > 6, A(r) always has a (single) positive real root for all values of a and therefore there exist 
black holes (which are always non-extremal) with arbitrarily large angular momentum per unit 
mass. These black holes were dubbed "ultraspinning" in ^46j and they will be the object of our 
study in this paper. 

In the limit of large angular momentum, these black holes can be characterised by two widely 
separated length scales on the horizon |46| . Let ^\\ denote the characteristic length scale in the 
directions parallel to the rotation plane, and ^±_ in the directions perpendicular to it. For a ^ r+, 
we have i\\ ~ a and ~ r+, that is, ^> Therefore, the MP black hole spreads out along 
the plane of rotation and it resembles a black membrane This observation has far reaching 
consequences because black branes are known to be unstable [47J. This led the authors of 06] 
to conjecture that rapidly rotating MP black holes should be unstable under Gregory-Laflamme 
type perturbations that preserve the symmetries of the {d — 4)-sphere. More precisely, choosing a 
cylindrical basis in polar coordinates (a = asvuO, (p) on the plane of the membrane, in the a — )• oo 
limit, the unstable modes should be of the form 



h^u-e^'JmiKa)e'"'^hf,,ir), (2.5) 

where k is the wavenumber along the direction parallel to the rotation. In this paper we will 
only address the axisymmetric case (m = 0), for which the radial profile is given by a cylindrical 
wave Jo(ko"). By extrapolating these observations to finite (but sufficiently large) o, one concludes 
that there should exist unstable modes that depend only on r and 6 and do not break any of the 
symmetries of the background black hole [46j . 

By considering the thermodynamics of singly-spinning MP black holes, the authors of [l6] 
observed that, for d > 6, these objects exhibit two markedly different behaviors depending on the 
value of a. For instance, consider the temperature of these black holes: 

+ — ]■ (2-6) 



I m 



^Notice that for the black ring of the angular momentum is bounded from below, but otherwise it can be 
arbitrarily large. 

^Ref. |46] showed that by taking a — )■ oo, — >■ oo and ^ — )■ while keeping f'^^ — rf^^ /a^ and a — a smO 



fixed, the singly-spinning MP black hole metric (2.11 near the rotation axis 9 = Q, reduces to the metric of a black 
membrane. In particular, the spatial directions along the worldvolume of the membrane correspond to the directions 
along the original plane of rotation in \2.1\ . 
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Figure 1: Phase diagram of entropy vs. angular momentum, at fixed mass, for MP black holes in 
d > 6 illustrating the conjecture of j46j (see also 140]): at sufficiently large spin the MP solution 
becomes unstable, and at the threshold of the instability a new branch of black holes with a central 
pinch appear (A). As the spin grows new branches of black holes with further axisymmetric pinches 
{B, C, . . . ) appear. We determine the points where the new branches appear, but it is not yet known 
in which directions they run. We also indicate that at the inflection point (0), where d'^S/dJ'^ = 0, 
there is a stationary perturbation that should not correspond to an instability nor a new branch 
but rather to a zero-mode that moves the solution along the curve of MP black holes, as discussed 
in subsection 12.21 



For fixed mass, r_|_ is a monotonically decreasing function of a. Therefore, it follows from (2.6) 
that, starting from a = 0, the temperature decreases as a increases, as in the Kerr black hole case. 
However, at 

the temperature reaches a minimum and then it starts growing like ~ r^^ as a increases. This is 
the kind of behavior of the temperature of a black membrane as the rotation increases. This point 
was proposed to give an order of magnitude for the appearance of the instability [l6] . Notice from 



(2.7) that the membrane- like behavior occurs for a > r+ in all dimensions (see the first column in 
Table [l]). The main result of this paper, which builds on |23] . will be to show that this picture is 
indeed correct. 

Furthermore, it is known ^48j that the threshold zero- mode (|A;| = fee, F = 0) of the GL in- 
stability gives rise to a new branch of static non-uniform strings. This led the authors of [l6] to 
conjecture that the zero- modes of the ultraspinning instability should give rise to a sequence of new 
branches of stationary and axisymmetric black holes that preserve all the M x U{1) x SO{d — 3) 
isometrics of the singly-spinning MP black hole, and that have ripples along the polar direction 6. 
This conjecture was further refined in HQ] , which argued that this structure of an infinite sequence 
of lumpy black holes is in fact needed in order to connect singly-spinning MP black holes to black 
rings, black Saturns and other black objects which include multiple concentric rings that should 
exist in d > 6 (see Fig{T]). 
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2.2 Black hole thermodynamics and the ultraspinning instability conjecture 

In this subsection, we will refine the arguments of |46j that we have reviewed above, and we will 
present the closely related conjecture made in Ref. [23]. The basic statement of this conjecture is 
that classical instabilities associated with stationary zero-modes can only appear in a regime which 
we call ultraspinning. We analyze this regime for the MP family of solutions. We will also argue 
that a stationary zero-mode of the black hole can correspond: (i) to a change in the parameters of 
the solution, if the zero-mode can be predicted by the condition of local thermodynamic stability, 
or (a) to the threshold of a classical instability of the black hole, which in our case will correspond 
to the prediction of [46j. 



2.2.1 The ultraspinning regime 

As discussed before, it was conjectured in Ref. [I6j that rapidly-rotating MP black holes with a 
single spin in d > 6 were unstable for Gregory-Laflamme-type modes. An order of magnitude 



estimate for the critical rotation, Eq. (2.7), was given by considering the thermodynamic behavior. 



The estimate for brane-like behavior is actually a zero-mode of the thermodynamic Hessian 

which we readily generalized to accomodate several spins. The condition of local thermodynamic 
stability is the positivity of this Hessian. It was shown in Ref. |24) that —Sa/3 possesses at least 
one negative eigenvalue for any asymptotically flat vacuum black hole (only the Smarr relation and 
the first law are required). Hence all such black holes are locally thermodynamically unstable. In 
the MP family, this is the only negative eigenvalue for small rotations. For > 6, as the rotation 
increases, the black hole acquires an additional thermodynamic instability, and a corresponding 
Gregory-Laflamme type instability of the associated black branes (this is a refinement of the Gubser- 
Mitra conjecture [5T]). 

Ref. |52j (see \22\ [23j for the rotating case) showed that a zero-mode of the thermodynamic 
Hessian is also a zero-mode of the action, i.e. an on-shell stationary perturbation of the black hole. 
However, this type of zero-mode consists of an infinitesimal change in the asymptotic charges of the 



black hole solution, within the MP family. In the singly-spinning case, Eq. (2.7) signals a degenerate 
point for which an infinitesimal change in the angular momentum at fixed mass does not change the 
temperature or the angular velocity of the MP black hole (indeed, our metric perturbations preserve 
these quantities, as we discuss later). This corresponds to the inflection point in Fig. [Tj Notice 
that we are using the MP equation of state S{M, Ji) to identify the zero-mode. Hence it does not 
have to be related to the bifurcation to a new family of solutions, and to the instability commonly 
associated with a bifurcation. Only zero-modes of the black hole which are not associated with the 



degeneracy of the Hessian (2.8) can represent bifurcations and instabilities. 



Ref. [23] and the present work confirm numerically the conjecture of Ref. [l6], showing that 



a zero-mode appears at (2.7), i.e. when the Hessian (2.8) is degenerate. More importantly, it 



is confirmed that additional zero-modes which are not thermodynamic in origin occur for higher 
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rotations. The latter zero-modes are thresholds for classical instabilities of the black hole, as 
explicitly verified for MP codimension-1 solutions in Ref. |24j . 

Zero-modes which are thermodynamic in origin can be identified as the zero-modes of the 
simpler reduced Hessian 



as shown in Ref. ^24j. In the MP family, the eigenvalues of Hij are all positive for small enough 
angular momenta, at fixed M. However, as some or all of the angular momenta are increased, 
some eigenvalues of Hij may become negative. We define the boundary of the region in which 
Hij is positive definite to be the ultraspinning surface. Following Ref. |23j . we shall say that a 
given black hole is ultraspinning if it lies outside the ultraspinning surface. Hence, as one crosses 
the ultraspinning surface, the black hole will develop a new thermodynamic instability, and the 
associated black branes will develop a new classical instability. This is in addition to the usual 
Gregory-Lafiamme instability already present at small angular momenta. 

Ref. [23] conjectured that classical instabilities whose threshold is a stationary and axisymmetric 
zero-mode occur only for rotations higher than a thermodynamic zero-mode, i.e. in the ultraspin- 
ning regime. We emphasize that the conjecture gives a necessary condition for this instability, not 
a sufficient one. Notice that by axisymmetric we mean that the Killing isometry (or ^id^j,. for 
more spins) is preserved. 

The intuition leading to the conjecture is that modes of lower symmetry are usually the most 
unstable ones. For instance, the original Gregory-Lafiamme instability occurs for the s-wave of 
the transversal black hole. An additional classical instability will arise after a critical value of 
the rotation, and it will correspond to a p-wave of the transversal black hole. As the rotation is 
increased, higher order waves may become unstable. Now, if we consider a black hole, instead of 
a black brane, the s-wave and the p-wave are associated with the asymptotic charges, mass and 
angular momenta. Therefore they are associated to purely thermodynamic instabilities. Higher 
order waves, on the other hand, may become classically unstable as the rotation is increased, 
starting with the d-wave. Notice that these waves do not affect the asymptotic charges}^ 

The fact that the thresholds of classical instabilities should be associated with bifurcations to 
different black hole families highlights the connection between stability and uniqueness. See e.g. 



2.2.2 Zero-modes of the Myers-Perry family 

Let us now examine the particular form of Hij in the case of general MP solutions, with N = 
[{d— l)/2j spins. We use the expressions in Ref. |27j . The horizon area, related to the entropy by 

^In the codimension-1 case (odd-dimensional equal spins) studied in Ref. [24l, where there is a precise harmonic 
decomposition of perturbations in terms of scalar harmonics of CP^ , it was explicitly shown that the first zero- mode 
changed only the angular momenta, breaking their equality, and that the next zero-modes (associated with classical 
instabilities/bifurcations) did not affect the asymptotic charges. 




(2.9) 



Ref. 00]. 
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S = A/4., is 



A 



Ad-2 



N 



\{{rl + a?) , 



(2.10) 



i=l 



where e = 0, 1 for odd and even d respectively. The surface gravity k, related to the temperature 
by T = k/2ti, and the angular velocities on the horizon are given by 



1 



N 

1=1 ^ ' 



2n 



(2.111 



The asymptotic charges, the mass M and the angular momenta Jj, which uniquely specify a 
solution, are given by 



M 



Ji 



167r ' d-2 

The reduced thermodynamic Hessian can be explicitly derived. 



a,M. 



(2.12) 



Hi 



{d-2)-K 
Mk 




K 



(r2+a2)2 (r2+a2)2 2r+ k 



(2.13) 



where there is no sum over and we have defined Cl"^ = X^j^^f- The matrix is positive definite 
in the static case, aj = Vi, and also in (non-extremal) d = 4,5. 
In the singly-spinning case, say a = ai 7^ 0, we have 



Hn 



Hi 



2{d-2){d-3)TTr+{rl + a^' 
M [{d - 3)rl + {d - 5)a^]^ 

id-2)7T 
MKrl 



[{d-3)rl - {d-5)a^] 



(2.14) 



for (i,i)/(l,l). 



There is a single zero- mode in d > 6 which occurs precisely for (2.7). The associated eigenvector is 
5ii, so that the angular momenta which vanish in the background solution are not excited by the 
perturbation, i.e. the zero-mode keeps the black hole within the singly-spinning MP family. The 
ultraspinning regime occurs for rotations higher than Eq. (2.7). Ref. |23] and the present work find 
classical instabilities whose thresholds occur in that regime only. 
In the equal spins case, whose analysis is also simple, we have 



H 



{d-2)-n r r 



Mk 



K 



2r^ 



ir'i + a2)2 



1 

2ru 



+ 



(2.15) 



where Qij = 1 Vi, j. An eigenvector Vi of Hij must then be an eigenvector of Qij, which leaves only 
two options: the eigenvector is such that Vi = Vyi, or is such that Vi = 0. In the former case, 
there can be no zero-mode, since this would require 



rl - a2 
(rl + a2)2 K. 



2r4 



{r\ + a2)2 2r 



1 iVr?2 
+ 



iV = 0. 



(2.16) 
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which can be simphfied to o? = —{d — 3)r^/(2 — e), and thus has no solution for a. Hence no 
instabihty occurs for modes that preserve the equahty between the spins, which is consistent with 
the results of [2H |69]. However, the — 1 eigenvalues associated with Vi = do change sign 
once at \a\ = r+j^ Therefore, an ultraspinning region exists even though the angular momenta are 
bounded by extremality, and indeed a classical instability was found in d = 9 [24|. 



2.3 The eigenvalue problem 

According to the preceding discussion, we are interested in modes that preserve the M x U{1) x 
SO{d — 3) symmetries of the background MP black hole (2.1), depending only on the radial and 
polar coordinates, r and 6. Thus we consider the following general ansatz for the perturbed metric 

V0 



ds^ 



A(r) 



dt + asin'^ee^'^ ^ 



2 sin^ 



(r^ + a^)d(j) — ae 



-Su) 



dt 



1 2 



+ e-'^o [dr + 6x sin 9 dOf + S(r, 9) e^^"' d9^ + cos^ 9 e''* dnl^_^^ , 



(2.17) 



where {5uq, 6ui,5fio, Sfii, 5uj, 5xi S^} are small quantities that describe our perturbations and they 
are functions of (r, 9) only. Unfortunately a decoupled master equation that governs the gravita- 
tional perturbations, analogous to the Teukolsky equation [33] in the case of the Kerr black hole. 



is not known for (2.1). Therefore, in this paper we will solve numerically the coupled partial dif- 
ferential equations (PDEs) that govern the class of perturbations we are interested in. Choosing 
the traceless-transverse (TT) gauge, 



h''^ = and Vh^^ = , 
the variation of the Ricci tensor in vacuum gives the following equations of motion: 

(Ai/i)^, = -VpV^V - 2 V/V = ' 



(2.18) 



(2.19) 



where Aj;, is the Lichnerowicz operator. Following [23] (see also [23]), we will consider a more 
general eigenvalue problem 

{AlH)^,^ = -klh^y . (2.20) 
This problem arises when one considers the stability of a uniform rotating black string. 



dsstring = gfj.u dx^dx" + dz^ , 

under perturbations of the form, 

'^'^string ^ ^'-'string ~\~ G ^ h^i,(^X^dx^dx 



(2.21) 



(2.22) 



^Notice that Ref. |24| used a difFerent radial variable, related to the variable presently used by — + a^, so 
that the ultraspinning surface is at \a\ = f+/V2- 



This ansatz is equivalent to the one presented in 
numerics (see subsection 3.51. 



but we found that this new form was better suited for the 
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where g^i/ is the metric of the singly-spinning MP black hole (2.1). The same problem arises 
when one considers the quadratic quantum corrections to the gravitational partition function in 
the saddle point approximation [67j (see [22] for the application to the Kerr-AdS black hole.) 



Our reason to consider (2.20) instead of trying to solve directly (2.19) is that Mathematica 
has very powerful built-in routines to solve generalized eigenvalue problems of this form. Thus, 
our approach will be to look for solutions of ( 2.20[ ) (which will generically have kc 0), and vary 
the rotation parameter a until we find a zero- mode, i.e., a mode with kc = 0. Therefore, the 
solutions with fee 7^ will correspond to new kinds of GL instabilities and inhomogeneous phases of 
ultraspinning black strings (see also [22]). On the other hand, the kc = modes will correspond to 
asymptotically flat vacuum black holes with deformed horizons of the kind conjectured in |46l I40j . 
as shown in Fig. [T] 

Finally, notice that the ansatz in (2.17) is the most general one that respects the isometrics 
of the background MP black hole and is preserved under diffeormorphisms that depend only on 
{r,6). Ultimately, this is necessary and sufficient to guarantee that (2.20) forms a closed system of 
equations 



3 Boundary conditions and gauge fixing 

In the following subsections we will discuss in detail the boundary conditions that we need to 
impose on the metric perturbations in order to solve (2.20). In the present situation, we have to 
specify boundary conditions at the horizon, r = r_|_, at asymptotic infinity, r — t- 00, at the rotation 
axis = 0, and at the equator 6 = 7r/2. In the following, we shall discuss the appropriate boundary 
conditions at the boundaries of the integration domain. In the final part of this section we show 
that the threshold stationary axismmetric zero-modes obeying these boundary conditions cannot 
be pure gauge. 



3.1 Boundary conditions at the horizon 

We shall demand regularity of the metric perturbations on the (future event) horizon by 
imposing regularity of the Euclideanised perturbed geometry on "H+J^ This approach allows us to 
discuss more straightforwardly the perturbations in the temperature and in the angular velocity. 
For r ~ r-(_, we can write A(r) = A'(r+)(r — r^) + 0[{r — r^)'^] with A'(r-i-) > oj^so that the 

^For the same reasoning, if we were considering time dependent non-axisymmetric perturbations (that preserve 
the transverse (d — 4)-sphere) we would also need to excite metric components of the type htr, hte, h^r, h^e, etc. We 
leave this very interesting problem for future work. 

^"This is equivalent to transforming the metric into regular coordinates on (e.g., ingoing Eddington-Finkelstein 
coordinates) and requiring that metric perturbations are finite on T-C^ in the new coordinate system. 

^^Recall that for d > 6 the singly-spinning MP black hole cannot be extremal and therefore A'(r+) > holds for 
all values of the rotation parameter a. The d = 5 case is special because there is a bound on the rotation whose 
saturation gives a nakedly singular solution. Therefore we shall only consider non-extremal solutions with A'(r+) > 0. 
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near horizon geometry of (|2.1|) reads 



A' (r+) (r — r+) 



dr' 



(3.1) 



+E {r+,9)de^ + 



sm 



a 

rl + a2 



S(r+,e) 

This suggests the introduction of a new azimuthal coordinate 

a 



rdt 



+ rl cos" 61 



H 



with 



t 



2tt. Wick-rotating the time coordinate, 

■ ^ ~ A'(r+) 



-IT . 



r 



with Tff 



2ttTh 

and defining a new radial coordinate p as 



47r (r^ + a^) 



and r ~ r + 



Th 



r+ + 



(3.2) 



(3.3) 



(3.4) 



the Euchdean near horizon (r ~ r+) geometry of the background solution (2.1) can be written in 
a manifestly regular form: 



ds^ f« S (r+, 6) [p'^dr'^ + dp"] + S (r+, 0) d^" + 



a")^ sin" I 



'" + cos" 6* . (3.5) 



S(r+,0) 

Indeed, at the axis of rotation {6 = 0) we have an explicitly regular S"" with no conical singularity. 
Moreover, the polar coordinate singularity at p = can be removed by a standard coordinate 
transformation into Cartesian coordinates {x,y). Note that a conical singularity at p = is 
avoided because we have chosen the period of the original Euclidean time coordinate f to be 
/3 = 1/Th. To sum up, regularity at the horizon of the background solution requires that we 
identify (r, (p) ~ (r, (/> + 2-ir) ~ (r + /3, — i Vinfi)- Furthermore, this procedure identifies Q-h with 
the angular velocity of the black hole and Th with its temperature. 

The boundary conditions for the metric perturbations can now be determined demanding that 
h^ydx^dx^ is a regular symmetric 2-tensor when expressed in coordinates where the background 
metric is regular. To do this, introduce manifestly regular 1-forms, 



= p^dr = X dy — ydx , = p dp = x dx + y dy . 

In terms of these 1-forms, the metric perturbation near the horizon reads 

'S(r+,^) A'(r+) + 2r+ (r^ + a") ' 



(3.6) 



h^y dx^ dx^ 



S(r+, 
-S(r+, 
-S(r+, 



5vq — 2a sin 6 



S2 (r+,^) A' (r+) 



6pi dO^ + 



^'(r+) p2 

2^2 ■ 2 - 

- a sm 



4i a(r^ + a") sin" g 5uj ~ (3 7) 
S(r+,^)A'(r+) ^^^^ 



(5i/i (i<^^ + cos" ^" 5$ dJ7^_4 



S(r+,( 
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Regularity of the perturbed geometry then requires that 



(5x, -5a; = 0(p2) , 5uo - Sfj^o = 0{p^) , dp6ni, dpSui, dp5^ = 0{p) . 



(3.8) 



near p = 0. 

It is important to note that we have imposed regularity of both, and separately, the background 
metric g^i^ and the perturbed metric hp_^. This implies that perturbations obeying (3.8) preserve 
the angular velocity and temperature of the background black hole 
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The regularity analysis of the boundary conditions is not complete without checking that the 



boundary conditions (3.8) are consistent both with the eigenvalue Lichnerowicz equations (2.19) and 



with the TT gauge conditions (2.18). We have explicitly confirmed this consistency. That is, the 



first term in the series expansion of the eigenvalue Lichnerowicz equations vanishes after we impose 



(3.8). On the other hand, we can use the TT gauge conditions to express e.g., {6uo,6i'i,6^} as 



functions of {5fio, 6pi, 5oj, 5x] and their first derivatives. Again, the first term of a series expansion 



of these TT gauge conditions is consistent with (3.8). 



3.2 Boundary conditions at the axis of rotation 

We follow the same strategy as in the previous subsection and we first discuss the regularity of the 



unperturbed background geometry (2.1) at the axis of rotation (0 = 0), where 9^ vanishes. In the 

nd ge 

A(r) 



region near = 0, the background geometry (2.1) reads 
A(r) 



ds 



J.2 _j_ q2 



dr + 2a 1 



J.2 _j_ g2 



dtdcp + p + a^) {de^ + e^dcp'^) + "^^T^dr^ + r'^dQ'^ 



A(r) 



(d-4) • 
(3.9) 

This metric can then be cast in manifestly regular form by changing to standard Cartesian coordi- 
nates (x, y) in the 6 — (p plane. 

To find the boundary conditions that the metric perturbations must satisfy at = we re- 
quire that hfj^^dx^dx'^ is a regular symmetric 2-tensor when expressed in coordinates where the 
background metric is regular. Introducing the manifestly regular 1-forms, 



E'^ = ede = xdx + ydy, 
the metric perturbation reads 



hpy dx^ dx'^ 



T^ovQdt + 

J.Z _j_ 

2a 



2 r"^ + 



A(r) 



6pQdr + 



Ar) 



xdy — ydx . 



SxE'^dr 



(3.10) 



+ ^^2^ [(r2 + a^ + A(r)) 6iv - (r^ + a^) M + A{r)6uo] E'^dt ^^'^^^ 



+ (r^ + a^) [5fiide^ + duiB^dcj)'^] + 6<^> dfl 



(d-4) • 



^^Notice that we could have chosen difFerent boundary conditions by requiring regularity of the full metric = 
+ Perturbations obeying such boundary conditions would generically change the temperature and angular 
velocity of the background geometry and we shall not consider them here, but would necessarily lead to a singular 
seen as a two-tensor on the background g^i^. 



15 



Regularity of the perturbed geometry then implies 

6vi - 6^1 = 0(6*^) , de6x, dgdio, dedfiQ, deduo, ded^ = 0{6) . 



(3.12) 



Again we have explicitly checked that the boundary conditions (3.12) are consistent both with the 



eigenvalue Lichnerowicz equations (2.19) and with the TT gauge conditions (2.18) 



3.3 Boundary conditions at the equator 



Introducing a new polar coordinate x = cos9, we find that the metric (3.13) near x = {9 = vr/2) 
is given by 



A(r) - 



2a 



P + - A(r)] dtdcj) + 



"2 + a^)^ - a'^A(r) 



(3.13) 



+ 



Afr) 



dr^ + dx + x dQ 



2 jri2 



{d-4) 



Once more, this geometry can be put in a manifestly regular form by changing to Cartesian coor- 
dinates. 

As in the previous subsections, we now demand that the metric perturbation h^ydx^^dx'^ is a 
regular symmetric 2-tensor when expressed in coordinates where the background metric is regular 
at a; = 0. Introducing the manifestly regular 1-forms, 



E"" = xdx, = x'^dn 



{d-4) 



(3.14) 



the metric perturbation near x = reads 



[A(r) - a2] 6uo + 2a'^5u 



dr + 



2a 



+ [(r^ + + A(r)) 6u) - (r^ + a^) 6i^i + A{r)5i^o\ dtdcp 



+ 



6fio dr 



A(r) A(r) \ x 

Regularity of the perturbed geometry for x — t- then requires that 

5x = 0{x) , (5$ - = 0{x^) , d^du, d^SfiQ, d^di^o, d^Svi = 0{x) 



(3.15) 



(3.16) 



Again we have explicitly checked that the boundary conditions (3.16) are consistent both with the 



eigenvalue Lichnerowicz equations (2.19) and with the TT gauge conditions (2.18). 



3.4 Boundary conditions at the asymptotic infinity 

At spatial infinity, r — )• oo, we will require that the perturbations preserve the asymptotic flatness 
of the spacetime. This means that they must decay strictly faster than the background asymptotic 
solution which, near spatial infinity, reduces to Minkowski spacetime 
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^■^More precisely, consider a spacetime {M.,g) which contains a spacelike hypersurface Mext diffeomorphic to R" \ 
B{R), where B{R) is a coordinate ball of radius R. The spacetime is said to be asymptotically fiat if the induced 
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In the asymptotic region, the eigenvalue Lichnerowicz equations (2.19) reduce simply to D/i 



-k^h^i,- The regular solutions (there are also irregular solutions that grow exponentially) of these 



equations decay as 



h 



1 



-kc r 



0, 



(3.19) 



for some constant a > that depends on the particular metric component and the number of 
spacetime dimensions d. Therefore, for kc ^ our perturbations decay exponentially to zero at the 
spatial infinity and hence asymptotic flatteness is guaranteed. 

Ultimately one would be interested in studying modes with kc = 0, since these correspond to 
exact perturbative solutions to Einstein vacuum equations in d spacetime dimensions. In this case, 
one would have to worry about the fall off of the metric perturbations. However, in our numerical 
method we will never be able to find modes for which kc = exactly, and therefore the fall off 
conditions at infinity for these modes are not an issue for us. 



3.5 Imposing the TT gauge conditions and the boundary conditions 



We have to solve the Lichnerowicz eigenvalue problem (2.20) for the seven metric perturbations de- 



scribed in (2.17), namely {6fio,6fj,i,5x,Su},6i'o, 51^1,5^}, subject to the TT gauge conditions (2.18). 
The gauge conditions allow us to eliminate three functions in terms of the other four and their first 
derivatives. In this paper we choose to express {6i^o, 61^1,6^} as functions of {6fio,5fj.i,6x,Sio} 
and their first derivatives. Notice that this choice differs from that in [23], where the independent 
variables were chosen to be {6fj,o,5fii,6x,S^}- The reason is that the present choice allowed us to 
obtain good numerical results in d = 5,6. As discussed in subsection 3.4, for kc ^ our (regular) 
perturbations vanish exponentially for r — )■ 00. However, in addition to the exponential decay, one 



also has a power law behaviour as in (3.19). This apparently irrelevant extra power law decay 
seems to make all the difference for the stability and/or accuracy of the numerical code in d = 5, 6. 
For d = 7, 8, 9 we obtain exactly the same results as in 



To summarize, we solve the gauge conditions (2.18) for {Suq, 6ui, 6^} in terms of {5/io, ^l^i, ^X^ 
and their first derivatives. Making these substitutions in the full set of the perturbation equations 



(2.20), we find that only four equations remain of second order in {6fiQ,6fj,i,6x,^i^}- Explicitly, 



metric h on Mcxt and the extrinsic curvature K satisfy 



5ij = Ok{r ") ; 



where r is the radial coordinate in R", and we write / — Okif^) if / satisfies 



(3.17) 



(3.18) 
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these equations are 



(Al/i)„ = - 
{AlHU = - 
{ALh)ee = - 

a{ALh)tt + 



-klKe , 
-k^hge , 

(r^ + a^) sin^ 9 



+ a^) sin^ 6 ^ ^ 



(3.20) 



ahtt + 



+ + sin^ ( 



ht<i, + 



(r2 + a^) sin^ 6 



and they constitute our final set of equations to be solved. 

A non-trivial consistency check of our procedure is to verify that the final equations (3.20) imply 



the remaining equations, which are of third order in the independent variables. We have verified 
that this is indeed the case. 



We will solve numerically the final eigenvalue problem (3.20) using spectral methods (see for 
instance |53|). In order to do so, we find it convenient to introduce new radial and polar coordinates 



y 



1 



!± 

r 



cos , 



(3.21) 



so that < y < 1 and < x < 1. The implementation of the method requires less computational 
power if all functions obey Dirichlet boundary conditions on all boundaries. Therefore, we redefine 
our independent functions according to 



ii{y,x) 
q3(.y,x) 



!± 

r 



x{l - x) 5fio{y,x) 
x{l - x) 6fj,i{y,x) 



q2{y,x) = r„^(l - x)6x{y,x) 
qi{y,x) = x6uj{y,x) , 



(3.22) 



so that the g^'s vanish at all boundaries. This guarantees that the boundary conditions discussed 



in the previous subsections, equations (3.8), (3.19), (3.16) and (3.12), are correctly imposed 



3.6 Zero-modes are not pure gauge 



In this subsection, we will argue that the TT gauge conditions (2.18) plus our boundary conditions 



(3.8), (3.19), (3.16) and (3.12), ensure that our modes are physical and not a gauge artifact. 



For kc > 0, the TT conditions completely fix the gauge since the action of the Lichnerowicz 
operator on a pure gauge mode is trivial, AiV(^^,^) = 0. However, TT perturbations with kc = 
can be pure gauge [17j and therefore it could be that the stationary perturbation fee — >• marking 
the onset of a new ultraspinning instability is not physical. In the rest of this subsection, we will 



show that there is no regular pure gauge perturbation that obeys our boundary conditions, (3.8) 



(3.19), (3.16) and (3.12). 



Under a gauge transformation with gauge parameter the metric perturbation transforms as 



(3.23) 
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The most general gauge parameter that preserves our ansatz (2.17) is of the form 



^ = (,r{r, x) dr + ix{r, x) dx , 



(3.24) 



where x = cos 6. We wih now prove that such a gauge parameter cannot generate a pure gauge 
metric perturbation that is regular on all boundaries. 

A TT gauge perturbation generated by must satisfy V^^'^ = and n^'^ = 0. If we introduce 
the antisymmetric tensor F^i, = V[^^,^] and we consider Ricci flat backgrounds, these conditions 
reduce to 

d^.{V^e)=0, d^{^gF^^)=Q. (3.25) 
Assuming that the solutions of these equations are separable, 



ir{r,x) = Rr{r)Xr{x) , Cxir,x) = R^{r)X^{x) 



the first equation in (|3.25|) reduces to 

= Ar'^-^ii,., dJx''-\l-x^)X., 



(3.26) 



(3.27) 



where A is a separation constant. Similarly, a combination of the r and x components of the second 
equation in (3.25) yields 



K r Rx , 



dAx''-\i - x')dxX, 



(3.28) 



where k is a second separation constant. Combining the equations (3.27) and (3.28) we further 
find that ^ ^ 

Rf = — R'x ) Xx = — X'j. . (3.29) 



We will take equations (3.27) and (3.29) as our independent equations. 
Plugging (3.29) into (3.27) gives 

dr (r'^-^AdrRx) = Kr'^-^Rx , dx fx^~^(l - x'^)dxX, 



The second of these equations has the solution 

fd-3-K d-3-K d-3 

Xr{x) = CiP 



7-d-K 7-d-K 7-d 



(3.30) 



, (3.31) 



where Ci^2 are integration constants and we have defined K = \J {d — 3)^ + 4k. For d > 6 (where 
we find the unstable modes) this solution diverges as C2X^~'^ at x = unless we set C2 = 0. In 
addition, a log(l— x) divergence at x = 1 can only be avoided if k satisfies the following quantization 
condition: 

K = 2^((i-3 + 2£), for ^ = 0,1,2,..., (3.32) 
which follows from the property r(— ^) = 00 for non-negative integer ^. 
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At this point we just have to solve the first equation in (3.30) with k given by (3.32). For our 
purposes it is sufficient to obtain the behavior of the solution near the horizon and near infinity. 
Near the horizon, r ~ r+, the most general solution behaves like 

„ , . Aai 1 



Rx{r) ^ao + ai log 1 



Rr(r) 



(3.33) 



2e{d - 3 + 2£) r - r+ ' 

where we used ( 3.29| ). However, equation (3.23) implies that such a diffeomorphism would generate 
metric perturbations of the form 

, _ ^Qi 1 , _ XaiPjx) 

'^'^ ^ 2£{d - 3 + 2£) (r - r+)2 ' ^ 2(r - r+) ' ^"^'"^^ 

which diverge at the horizon. Therefore, our regularity requirements at the horizon force us to set 
oi = 0. Similarly, near infinity a Frobenius analysis gives the asymptotic solution 

i?,(r)«6or-('-'+'') + 6ir2^ (3.35) 



The term in this equation proportional to 6i generates, through (3.23), a metric perturbation 
of the form hrr oc bir'^^^~^^ near infinity. This perturbation must decay strictly faster than the 
unperturbed background solution (which at infinity is Minkowski), which implies that we have to 
impose bi = for any £ >0. 

Summarizing, the gauge parameter (3.24) that could potentially generate a metric perturbation 

(3.36) 



that is regular both at the horizon and at infinity must have the asymptotic behavior 
We can now complete our proof. Notice that 



< 



dr/-^A(r) [drRxir)Y 



A{r)R,{r)drRx{r] 



drdr 



/>oo 

-2i{d-3 + 2£) /"^i?x(r)^ < 0, 



A{r)drRAr) R^r) 



(3.37) 



where we used (3.36) and (3.28). But these relations can be satisfied only for Rx{r) = 0, which in 
turn implies that Rr{r) = and hence = 0. 

Therefore, we have proved that there is no regular gauge parameter (3.24) that could potentially 
generate the metric perturbations that we consider, eq. (2.17). Thus, we conclude that our regular 
zero-mode perturbations cannot be pure gauge. 



4 Results 

In this section, we present our results for the spectrum of negative modes of the Lichnerowicz 
operator. The actual spectrum in d = 5, 6, 7 is displayed in Fig. [2j for the other values of d up to 
d = 11 the results are qualitatively similar to d = 6, 7 and we expect the same to be true for all 
values of d with d > 11. Following [23], we plot the dimensionless negative eigenvalue —k'^r'^ as 
a function of the (dimensionless) rotation parameter o/r^- In Figs, [s] and [4] we show the actual 
metric perturbations {dpLQ^ 5/Ui, 5x, 5w} for d = 7, and in Fig. |6]we display the embeddings of the 
unperturbed and the perturbed horizons for the same number of spacetime dimensions. 
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Figure 2: Negative modes of the singly-spinning MP black hole m d = 5 (left), d = 6 (centre), and 
d = 7 (right). 



d 


ia/rm)\e=i 


ia/rm)\e=2 


ia/rm)\i=3 


{a/rm)\i=4 


6 


1.097 


1.572 


1.849 


2.036 


7 


1.075 


1.714 


2.141 


2.487 


8 


1.061 


1.770 


2.275 


2.725 


9 


1.051 


1.792 


2.337 


2.807 


10 


1.042 


1.795 


2.361 


2.855 


11 


1.035 


1.798 


2.373 


2.879 



Table 1: Values of the rotation a/r^ for the first four harmonics of stationary perturbation modes 
{kc = 0). The estimated numerical error is ±3 x 10"'^ in d = 6, 7 and ±5 x 10^'^ in d = 8, 9, 10, 11. 



4.1 Results for d = 5 

Ref. [21] proved that MP black holes (not necessarily singly-spinning) are always locally thermody- 
namically unstable. Therefore, it follows from the arguments in subsection |2.2| that the spectrum 
of the Lichnerowicz operator should admit at least one negative eigenvalue. Since in d = 5 the 



Hessian (2.14) is positive definite, the expectation is that the Lichnerowicz operator should admit 
one and only one negative mode. Our results for the spectrum of the Lichnerowicz operator are 
depicted in Fig. [2] (left) and indeed they confirm this expectation. As the rotation is increased, kc 
increases, in agreement with the results of \Q9\ [22l [23] . If we interpret this result from the black 
string perspective, then more modes < kc are GL-unstable, which is expected since the centrifu- 
gal force should make the string "more unstable". The negative mode diverges as kc oc l/{rm — a), 
in the singular limit a — >• r^- 



4.2 Results for 6 < d < 11 

The situation in d > 6 is more interesting. The spectrum of the Lichnerowicz operator in c? = 6, 7 
is displayed in Fig. [2] (center) and (right). We label the different branches that intersect kc = 
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at finite a/rm by successive integers I = 1,2,3,..., and we refer to the corresponding modes as 
harmonics, although the equations we are solving do not seem to separate. The values of a/rm at 
which the first few branches intersect fee = are summarised in Table [T] We note that for each 
branch the corresponding integer £ coincides with the number of zeros that the metric perturbations 
Siio{x,y) and 5jii{x,y) have on the horizon y = (see Fig. [s}) 

First, we notice that the Lichnerowicz operator has a negative eigenvalue for all values of 



a/rm- This is in agreement with the thermodynamic argument in subsection 2.2, However, for 
d > 6, singly-spinning MP black holes admit an ultraspinning regime. The ultraspinning surface, 
which determines the boundary of this region (in parameter space), is given by the onset of the 



membrane-like behavior (see eq. (2.7)): 

aV'' d-3 fd-^V'-"^'^ 

(4.1) 



^-/mem 2(d - 4) " 5 



According to the thermodynamic argument in subsection 2.2, precisely at this value of a/rm the 
Lichnerowicz operator should develop a new negative eigenvalue and we find that this is indeed the 
case. In the first column of Table [T| we display our numerical results for the critical values of a/rm 
for the £ = 1 mode in various dimensions. We find that the numerical results agree very well with 



the values oi a/rm computed from (4.1 ). We emphasise that the £ = 1 mode should correspond to a 
variation of the parameters within the MP family of solutions such that it preserves the temperature 
and angular velocity of the background (this mode is in correspondence with point in Fig. [T]). 
Therefore this mode should not correspond to an instability of the black hole. However, it should 
give rise to a new type of classical instability of the associated black string, in which the horizon is 
deformed along both the direction of the string and the polar direction of the sphere. 

For a/vm > (a/?"m)mem, 1-6., in the ultraspinning regime, there appears an infinite sequence of 
new negative modes. The values of a/rm at which the i = 2,3,4 branches intersect kc = are 
displayed in Table [l} These modes do not admit a thermodynamic interpretation and therefore 
they should correspond to new perturbative black holes with deformed horizons (see points A, B, C 
in Fig. [T]). In particular, the £ = 2 mode should signal the onset of the ultraspinning instability 
conjectured in |46j. For d = 7 the actual metric perturbations 6fio, ^Aii, Sx, for the £ = 1,2 
modes are displayed in Figs. [3]j4j Figures 3(a)[|3(b~ 3(c) and |3(d)l show that for each of the ^ = 1, 2 



modes, the number of zeros that Sfj,o{x, y) and dfj,i{x, y) have on the horizon [y = 0) coincides with 
the integer £. We have checked that the same is also true for the higher harmonics with £ > 2. 



According to the discussion in subsection 2.1, an important prediction of [46J is that, in the 



a — 7- oo limit and in a region close to the axis ~ 0, the threshold axisymmetric modes at the 
horizon should be well-approximated by a Bessel function Jo(ko"). In order to check if, for large 
a and close to the horizon and the axis, our perturbations reduce to a Bessel function (at least 
qualitatively), we have fitted our numerical results in d = 7 for A(r) hrr with a/rm = 4.38 and 
^c'^m = 0.114j^ with a Bessel function: 

aJo(Kcr), with a = -1.868, k = 5.564, (4.2) 

"For these results y = 0.0034. 



22 




Figure 3: Functions 5fj,o{x,y) and 5^i{x,y) for the i = 1 (Figs. 3(a) and 3(c)) and £ = 2 (Figs. 
|3(b)| and |3(d)[ ) modes respectively. The number of zeros at y = (the horizon) coincides with the 
integer £. 



and a = a sin 9. The results are depicted in Fig. [5] and they show that the agreement is quite 
remarkable. We should emphasise though that the argument of |46j only applies in the strict limit 
mentioned above. For our data a/r^ is relatively large (compared to the onset of the ultraspinning 
instability) but nevertheless finite, and therefore only a qualitative agreement with the prediction 
of |46| should be expected. This is what Fig. [5] shows. 

To confirm our interpretation and visualise the effect of these perturbations on the horizon of the 
background MP black hole, in Fig. [6] we compare the embedding diagrams of the unperturbed and 
the perturbed horizons in d = 7 (see Appendix [A| for the technical details of the embedding diagram 
construction). For other values of d the picture is qualitatively similar. From the embeddings we 
conclude that the effect of the i = 2 modes is to create a pinch centered on the axis of rotation; 
the i = 3 modes create a pinch at a finite latitude and the 1 = 4 modes create two pinches, one 
centered on the rotation axis and the other at a finite latitude. These are precisely the kind of 
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(a) 



(b) 





0.0 0.0 



(c) 



(d) 



Figure 4: Functions Sxix,y) and 5uj{x,y) for the 
and |4(d)[ ) modes respectively. 



1 (Figs. [4(^0] and [4(H)]) and ^ = 2 (Figs. [4(b)] 



deformations depicted in Fig. [T] 

Finally we notice that the value of (a/rm)crit for every i > 1 increases with the number of 
spacetime dimensions d. It would be interesting to explore if the perturbative approach of [17j can 
capture the dynamics of these instabilities in the a — >• oo limit. 



5 Discussion 

In this paper we have studied in detail the properties of the onset of the ultraspinning instability 
in singly-spinning MP black holes in d > 6. This instability is captured by a class of perturbations 
that preserve the M x C/(l) x SO{d — 3) symmetries of the background, as well as the angular 
velocity and temperature of the original MP black hole. These perturbations (in the TT gauge) 
must satisfy the Lichnerowicz system of equations (2.19). Our strategy however was to solve 



the more general eingenvalue problem (2.20) because: i) there are powerful numerical routines to 
solve generalized eigenvalue problems of this kind; and ii) this also allows us to study Gregory- 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

0"/rm 



Figure 5: Comparison between our A{r)hrr for a/rm = 4.38 and fe^r^ = 0.114 in d = 7 (red 
dots) and the fitting function aJQ{Ka) (dashed hne). The red line corresponds the interpolating 
polynomial to our data and should only serve to guide the eye. To obtain the plots, the minimum 
value of X was taken to be Xmin = 0.69804, since the comparison should be relevant near the rotation 
axis at X = 1. 



Laflamme— like instabilities of the associated black string. In short, varying the rotation parameter 
a of the background black hole (2.1 ) we searched for the negative modes of the problem (2.20), i.e., 
for solutions that have generically kc 0. We found several families of negative modes that exhibit 
an underlying harmonic structure, although the equations that we solve do not seem to separate. 
We used this to suggestively label the several branches by an integer i. This integer coincides with 
the number of zeros of the metric perturbations on the horizon. The family i = always has 
fee 7^ 0, and coincides with the well-known negative mode of the Schwarzchild-Tangherlini solution 
when the rotation parameter a vanishes [67, 70] . This family therefore describes the critical mode 
I A; I = kc of the original Gregory-Laflamme instability of the black string |47j when a = 0, and our 
results show how it evolves as the rotation increases: the value of the threshold wavenumber kc 
increases with the rotation. 

The branches with £ > 1 are more interesting since they intersect the kc = axis at a critical 
a > (see Table [T]). Thus they describe not only new types of Gregory-Laflamme instabilities of 
rotating black strings, but they also represent, for i >2, true instabilities of the MP black hole. For 
black strings, the onset of these instabilities is conjectured to signal a bifurcation to new branches 
of non-uniform black strings in which the horizon is deformed along both the direction of the string 
and the polar direction of the transverse sphere. At a given rotation, the i = 1 mode has the 
shortest wavelength k~^ and hence it should dominate the instability of the black string. 
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From now on, we focus our attention only on the consequences of our findings for the stabiUty 
of the MP black holes. The critical rotation where the i = 1 mode appears can be predicted by 
the thermodynamic argument of \23\ I24j (see Section 2.2), which can be seen as a refinement of the 
Gubser-Mitra conjecture. Accordingly, this mode does not correspond to a true instability of the 
MP black hole. Instead, we have interpreted it as a thermodynamic mode that corresponds to a 
variation of the parameters within the MP family of solutions such that it preserves the temperature 
and angular velocity of the background (Section [4| and point in Fig. [T]). 

The modes with £ > 2 describe the onset of true ultraspinning instabilities of the MP black 
hole. In section |4j we have studied the deformations that these modes produce on the horizon of the 
black hole for d = 7 (just for concreteness) and found that they give rise to the kind of deformations 
predicted in [401 . In addition, we have shown that, for large a and near ^ = 0, our numerical results 
are well- approximated by a Bessel function, which is in agreement with the heuristic arguments 
of [l6]. Altogether, these results provide solid evidence in favor of the ultraspinning instability 
conjectured in \A6\ I40j . The thresholds of the ultraspinning instabilities are expected to signal 
bifurcations to new branches of axisymmetric solutions with pinched horizons in a phase diagram 
of stationary solutions (see points A,B,C,... in Fig. [T]). Although these pinched black holes have 
the same isometrics as the original MP black hole, their spherical-topology horizons are distorted by 
ripples along the polar direction. When continued in the full non-linear regime, these new branches 
of solutions are conjectured to connect to the black ring and black Saturn families [40j. Although 
we have identified the critical values of the rotation where these new branches of pinched black 
holes appear, unfortunately with our methods we cannot determine if they will have larger or lower 
entropy in the phase diagram S vs. J at fixed total mass. This would require going beyond linear 
order in perturbation theory. 

We have limited our study to stationary perturbations and therefore we cannot claim to have 
found an unstable mode, i.e., a linear perturbation (satisfying our boundary conditions) that grows 
exponentially with time; we have only found the stationary zero-modes that signal the onset of 
the instability. Including time-dependence is not conceptually difficult but it is only technically 
harder. As pointed out in footnote [9], consistency of (2.20) requires turning on extra components 
of the perturbed metric. This problem is however of fundamental interest since its solution would 
provide a definite proof of the ultraspinning instability together with information on its timescale. 
The analogous stability problem, including time-dependence, was studied in [21] in the context of 
odd-dimensional MP black holes with equal angular momenta in all rotation planes. There, the 
analogue of the £ = 1 zero-mode is also present, and the time-dependent analysis confirmed the 
absence of a black hole instability in this sector. We take this result as good evidence in support 
of a similar interpretation in the singly-spinning case. 

The d = 5 case is special. As |46j already pointed out, the d = 5 singly-spinning MP black hole 
does not have an ultraspinning regime (in the sense of subsection 2.2) and therefore a priori there 
is no argument that suggests the existence of an instability within the class of perturbations that 
preserve the isometries of the background. In this paper we have confirmed this picture since in 
d = 5 we only found one negative mode (the £ = branch) which has kc ^ for all values of a such 
that a < Vm- The d = 4 case has already been discussed in [22] . 

An interesting open question to be addressed is whether the ultraspinning instability is also 
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present in MP black holes in AdS. This background might introduce new features and physics 
that deserve a study. In particular, AdS might exhibit new black hole phases with an interesting 
interpretation in the holographic dual field theory. In fact, pinched plasma balls [Mj, new kinds 
of deformed plasma tubes [62] and rotating plasma ball instabilities |65] have been found in the 
context of Scherk-Schwarz-AdS. The spectral methods that we have employed can also be applied 
in AdS spacetimes (see [22J for an application in d = 4) . 

The existence of an ultraspinning regime is not unique to MP black holes. For instance, the five- 
dimensional black rings also admit a certain ultraspinning regime. Namely, in the limit of arbitrarily 
large angular momentum the black ring resembles a boosted black string [BU], which is known to 
suffer from GL instabilities |54j . Moreover, the arguments in |28l [5^ [55] suggest that black rings 
should suffer from this GL instability for relatively low values of the angular momentum. Although 
these GL modes would break the symmetries of the background, it is interesting to ask whether 
there are zero-modes which can be captured by our methods. In fact, black rings are expected to 
suffer from other types of instabilities. For instance, fat rings are expected to be unstable under 
changes of their radius |55j . and doubly-spinning black rings are conjectured to suffer from similar 
instabilities as well as from superradiant instabilities associated to the rotation along 5^ |56j . 

Rotating asymptotically flat black holes can also suffer from instabilities which break the ax- 
isymmetry. Recently, \19\ [20] performed a remarkable fully non-linear numerical evolution of a 
perturbed singly-spinning MP solution in d > 5 and found that these black holes are dynami- 
cally unstable against non-axisymmetric bar-mode perturbations. The endpoint of this instability 
(at least in the regime of rotations explored in [20]) seems to be another MP black hole with a 
smaller angular momentum. Notice however that, because the axisymetry is broken, the threshold 
of the instability is not expected to be stationary, and thus there will be no bifurcation to a new 
stationary family of black holes. The existence of such an instability had been conjectured also 
by Emparan and Myers in [l6] using a thermodynamic argument of a different type to the one 



explored in the current manuscript This bar-mode instability is certainly different in nature 
from the axisymmetric instability, which is not present in the d = 5 case. Moreover, the bar-mode 
instability becomes active at slightly lower rotations: compare the critical values of a/r^ for the 
onset of the bar instability in Table 1 of [2D] with the critical values of a/r^ for the onset of the 
ultraspinning instability in our Table [T] It would be certainly very interesting to study also the 
non-linear time evolution of the axisymmetric ultraspinning instability. 

In addition, in the presence of certain matter or a cosmological constant, rotating or charged 
black holes (including the Kerr solution) can be unstable due to the superradiant phenomena 
[501 [571 [58l [59] . Consider a black hole with angular velocity (or with chemical potential fJ-n)- 
If the frequency oj and the angular momentum m (or charge e) of a scattering wave satisfy the 
relation u < m^^n (or w < efin), the scattering is superradiant: the wave extracts rotational 
energy (or electromagnetic energy) from the background and gets amplified. The presence of an 
effective reflective potential barrier, which may be due to the mass in a massive scalar field case 
or due to the AdS gravitational potential, combined with superradiance can then lead to multiple 



^^The original argument of [46] for the bar-mode instability is the following. The entropy of a highly rotating black 
hole can be smaller than that of two boosted Schwarzchild-Tangherlini black holes in orbital motion. Therefore it 
should be entropically favored to the system to develop a non-axisymmetric configuration. 
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wave amplification and reflection that renders the system unstable. This can occur for scalar, 
electromagnetic or gravitational (but not fermionic) waves. Such an unstable system is often 
coined as a "black hole bomb" [57J . It would be very interesting to determine the endpoint of this 
superradiant instability \50\ [58] and there is active ongoing research in this direction [60j. The 
spectral methods might also be useful to tackle this problem. 

Returning back to Smarr's expectation of an instability in rotating black holes based on an 
analogy between rotating black holes and fluid droplets |44j . it is interesting that it is realised in 
higher dimensions with the ultraspinning instability. That this happens is no longer seen as a mere 
analogy but a consequence of the gauge/gravity correspondence. Black holes are thermodynamic 
objects and as such one should expect that, in a long wavelength regime (compared to the energy 
scale set by the temperature), they should have an effective hydrodynamic description. The first 
attempt to materialize this idea was the membrane paradigm |Blj. More recently, the AdS/CFT 
duality motivated further research in this direction which culminated with the formulation of the 
fluid/gravity correspondence [12], which provides a precise hydrodynamic description of asymptot- 
ically large AdS black holes. In addition, we should mention the blackfold approach of |17| [TH] 
which also provides a hydrodynamic description of asymptotically flat black holes in the regime 
where this approximation applies. In both latter cases, the formal connection between gravity and 
fluid dynamics is established through a derivative expansion of the Einstein equations. 

In this context, it is therefore with no surprise that many of the aforementioned black hole 
instabilities have an hydrodynamic description. The Gregory-Laflamme instability |47j of black 
strings is in correspondence with the Rayleigh-Plateau instability in fluid tubes j62j and with 
damped unstable sound wave oscillations [63]. The subject of our study, the ultraspinning instability 
and the associated new phases of pinched stationary black holes, also possesses a fluid description 
as ultraspinning pinched plasma balls [|64j . Finally, the bar mode non-axisymmetric instability [19] 
was also conjectured to exist in black holes due to its presence in rotating fluids [65] . 

Acknowledgements 

It is a pleasure to acknowledge the stimulating discussions with our collaborators Roberto Emparan 
and Harvey Reall in the ultraspinning project. OJCD acknowledges financial support provided 
by the European Community through the Intra-European Marie Curie contract PIEF-GA-2008- 
220197. PF is supported by an STFC rolling grant. RM and JES acknowledge support from the 
Fundagao para a Ciencia e Tecnologia (FCT-Portugal) through the grants SFRH/BD/22211/2005 
(RM) and SFRH/BD/22058/2005 (JES). This work was partially funded by FCT-Portugal through 
projects PTDC/FIS/099293/2008, CERN/FP/83508/2008 and CERN/FP/109306/2009. 

A Embedding diagrams 

In this appendix, we summarise the main ingredients of the embeddings presented in Fig. [6j 
Since the kind of perturbations we are considering preserve the transverse {d — 4)-sphere, we will 
henceforth suppress it and concentrate on the non-trivial four-dimensional part of the metric. In 
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order to be able to visualise the geometry of the horizon for all values of the rotation parameter 
we will adopt the embedding proposed in [7T] in the context of the Kerr-Newman black hole. 

Recall that for both the perturbed and the unperturbed geometries, the spatial cross-sections 
of the future event horizon are hypersurfaces of constant t and r = r^. Here r-(_ denotes 
the location of the event horizon of the background MP black hole. The induced metric on this 
hypersurface (suppressing the transverse {d — 4)-sphere) can be written as 



2 jj,2 



with a 



, and V{a) = 



dx\ 



(A.l) 



where x = cos 6. This two-dimensional surface of revolution can be embedded into Euclidean 
four-dimensional space via the map (a, 0) i— )• {X, Y, Z, R) defined by 

{X,Y,Z) = ^{F{<P),G{^),H{^)), R = R{a), and F{4>f + G{<t>f + H{(pf = I , (A.2) 



and the induced metric for the associated 2-surface is 



dsh = dX^ + dV^ + dZ^ + dR^ 



-2 



dR\- 
d^J 



da^ + a^dd)^ 



(A.3) 



where we ask the reader to see [7T] for a detailed understanding of the second equality. 
Introducing the scale parameter ij and the distortion parameter /3, 



T] 



it follows that 

9xxl.^ = rff{x) , gcf>^\,^.^ = ri^f(xy^ , with f{x) = ^ _^2(i"_^2) 
Then, from ds'^4 = ds\^, we get 

Reality of this expression requires that $0 ^ \ f'{x)\/2 for < x < 1 which is satisfied if 

^o>\ max . 
2 xe[0,l\ 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



In the embedding diagrams shown in Fig. |6j the Cartesian coordinate Z was chosen to be given 
by the R{x) that solves (A.6). Similarly, X was chosen to be given by the function /(x)^/^ defined 
in 
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Figure 6: Embedding diagrams (Figs 6(a)[6(b) 6(c) ) at {a/rm)cnt of tlie d = 7 black hole horizon, 
unperturbed (blue), and perturbed (red) with the i = 2,3,4 harmonics (in this order). The 
embedding Cartesian coordinates Z and X lie along the rotation axis 9 = and the rotation 
plane ^ = 7r/2 respectively. We also show the logarithmic difference between the embeddings of 
the perturbed (^^=2,3,4) and unperturbed (Zq) horizons. The number of spikes corresponds to the 
number of crossings between the two embeddings. Each red + sign indicates that the perturbation 
bulges out relatively to the background, and the blue — signs indicates the opposite situation. 
Figures |6(d) 6(e) and 6(f)| provide a three-dimensional illustration of the spatial sections of the 
horizon, with the (d — 4)-sphere suppressed at every point, with the same color code. 
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